Complexity of waves in nonlinear disordered media 
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The statistical properties of the phases of several modes nonlinearly coupled in a random system 
are investigated by means of a Hamiltonian model with disordered couplings. The regime in which 
the modes have a stationary distribution of their energies and the phases are coupled is studied for 
arbitrary degrees of randomness and energy. The complexity versus temperature and strength of 
nonlinearity is calculated. A phase diagram is derived in terms of the stored energy and amount 
of disorder. Implications in random lasing, nonlinear wave propagation and finite temperature 
Bose-Einstein condensation are discussed. 



The interplay between disorder and nonlinearity in 
wave-propagation is a technically challenging process. 
Such a problem arises in several different frameworks 
in modern physics, as nonlinear optical propagation and 
laser emission in random systems, Bose-Einstein conden- 
sation (BEC) and Anderson localization (as, e.g., in Refs. 
[Il49ll9l4l9l|). Related topics are the super-continuum gen- 
eration and condensation processesi^— 

When disorder has a leading role, nonlinear processes 
can be largely hampered as due to the fact that waves 
rapidly diffuse in the system. Conversely, if the struc- 
tural disorder is perturbative, its effect on nonlinear evo- 
lution is typically marginal, leading to some additional 
linear or nonlinear scattering-losses, but not radically af- 
fecting the qualitative nonlinear regime expected in the 
absence of disorder. When disorder and nonlinearity play 
on the same ground, one can envisage novel and fascinat- 
ing physical phenomena; however, the technical analysis 
is rather difficult, as the problem cannot be attacked by 
perturbational expansions. 

Physically, disorder and nonlinearity compete in those 
regimes when wave scattering affects the degree of local- 
ization, eventually inducing it (as in the Anderson lo- 
calization), and nonlinearity couples the modes in the 
system. These may in general exhibit a distribution of 
localization lengths (determined by the amount of disor- 
der) and a strength of the interaction depending on the 
amount of energy coupled in the system. 

Our interest here is to provide a general theoretical 
framework, whose result is the prediction of specific tran- 
sitions from incoherent to coherent regimes, which are 
specifically due to the disorder and display a glassy char- 
acter, associated with a large number of degenerate states 
present in the system. 

We adopt a statistical mechanics perspective to the 
problem, which allows to derive very general conclusions, 
not depending on the specific problem, and our focus is 
on the case in which many modes are excited. This im- 
plies that energy is distributed among many excitations 
in an initial stage of the dynamics. The overall coherence 
(i.e., the statistical properties of the overall wave) will be 
determined by the phase-relations between the involved 
modes. Here we show that there exist collective disor- 
dered regimes, where coherence is dictated by the fact 



that the system is trapped in one of many energetically 
equivalent states, as described below. 

Representing mode phases by means of continuous pla- 
nar XY-like spins and applying a statistical mechanic ap- 
proach we can identify different thermodynamic phases. 
For negligible nonlinearity, all the modes will oscil- 
late independently in a continuous wave noisy regime 
("paramagnetic" -like phase). For a strong interaction 
and a suitable sign of the nonlinear coefficients, all 
the modes will oscillate coherently ("ferromagnetic" -like 
regime). This corresponds, for example, to standard pas- 
sively mode-locked laser systems 26 , which we found to 
take place even in the presence of a certain amount of dis- 
order. In intermediate regimes, the tendency to oscillate 
synchronously will be frustrated by disorder, resulting in 
a glassy regime. 

These three regimes are identified by a set of order pa- 
rameters (up to 10 for the most complicated phase, as 
detailed below), which can be cast into two classes: the 
"magnetizations" to, and the "overlaps" q. As the sys- 
tem is in the paramagnetic-like phase all to and q vanish; 
in the ferromagnetic regime they are both different from 
zero; while in the glassy phase to = and (at least some 
of) the overlap parameters are different from zero. 

The paramagnetic and ferromagnetic phases may be 
present even in the absence of disorder; conversely, a 
necessary condition to find a glassy phase is frustration 
(disorder induced in our case, see Sec. IIII II) . The glassy 
phase is characterized by the occurrence of a rugged - 
complex - landscape for the Gibbs free energy functional 
in the mode phases space: a huge number of minima 
are present, corresponding to a multitude of stable and 
metastable states in the system, separated by barriers of 
various heights and clustered in basins. This is a result 
of the competition between disorder and nonlinearity. 

The existence of a not-vanishing complexity (which 
measures the number of energetically equivalent states) 
for the possible distributions of mode-phases is the basic 
ingredient for explaining a variety of novel phenomena 
like speckle pattern fluctuations and spectral statistics 
for disordered, or weakly disordered nonlinear, systems, 
ergodicity breaking, glassy transitions of light or BEC, 
and ultimately the onset of a coherent regime in a ran- 
dom nonlinear system. 
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Our work extends previously reported results, cf. Ref. 
[27| and it includes an arbitrary degree of disorder and 
the discussion of its application to nonlinear Schroedinger 
models, relevant, e.g., for BEC, spatial nonlinear optics 
and supercontinuum generation. 

The paper is organized as follows: in Sec. U we intro- 
duce the model and we discuss some of its possible fields 
of application, namely random lasers, Bose-Einstein con- 
densates, optical propagation; in Sec. |TT]we discuss the 
effect of disorder in the coupling of light modes and the 
new expected phenomena; we dedicate Sec. IHII to an ex- 
tremely basic introduction to the statistical mechanics of 
systems with quenched disorder, to the replica method, 
and to the definition of complexity; in Sec. IIVI we study 
the model within the replica approach, details of the com- 
putation are reported in App. [A] in Sec. |V] we discuss 
the presence of excited metastable states and we com- 
pute the complexity functional; in Sec. IVII we show the 
phase diagrams of our model and discuss the properties 
of its thermodynamic phases; eventually, in Sec. IVIII we 
draw our conclusions. 



I. THE LEADING MODEL 

Here we review some of the disordered systems where 
a relevant non linear interaction may arise and our model 
applies. The basic Hamiltonian of N adimensional angu- 
lar variables <j) € [0 : 2it] is given by 



localization and not-negligible nonlinear effects are im- 
portant; a few of them are detailed in the following. 

As a thermodynamic approach is adopted, one can 
argue if the statistical mechanic techniques also apply 
in those systems where the definition of a temperature 
is not straightforward, as, specifically, nonlinear optical 
wave propagation in disordered media. This particular 
problem can, then, be treated as for constraint satisfac- 
tion problems in computer science^r— where - at the 
end of the calculation - the limit of zero-temperature is 
taken and it is shown that a transition is expected as the 
number of constraints grows. 



A. Random active and passive electromagnetic 
cavities 

We start from the electromagnetic energy inside a di- 
electric cavity (due to the generality of the considered 
model similar examples can be found in a variety of sys- 
tems): 

£em= /E(r)-D(r) dV (2) 

The displacement vector is written in terms of a position 
dependent relative dielectric constant £r(r): 

D(r) = e er(r)E(r) + e PAr L (r) (3) 
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where i = {ii,i2,i$ii} and J\ are random independent 
identically distributed interaction variables. Formally, 
the couplings can vary from short- to long-range, de- 
pending on the structure of the four-index interaction 
tensor J\. If we choose Ji =/= for any distinct quadru- 
ple ix, ■ ■ • ,i4, independently of the geometric position, we 
can build a mean-field theory in which the system is fully 
connected. In this case the average Ji and the variance 
of its distribution must scale as 1/N 3 to guarantee ther- 
modynamic convergence of (free) energy density. The in- 
teraction can, otherwise, be bond-diluted with arbitrary 
degree of diluteness, adopting a sparse tensor whose non- 
zero elements do not scale with the number of modes, 

As we show in the following, the Hamiltonian, Eq. (TT]) 
is derived in different contexts, and the various parame- 
ters may have different interpretation. In this manuscript 
we want to derive general properties that are expected 
assuming a simple, yet reasonable, Gaussian distribu- 
tion for the random coupling coefficients, with a non- 
vanishing mean value. Varying the ratio between the 
standard deviation and the mean value we control a dif- 
ferent degree of disorder. Hence, these results applies 
to the various cases in which random wave propagation, 



with Pnl the nonlinear polarization. In absence of the 
latter, for a closed cavity, the field can be expanded in 
terms of the modes of the system. In the presence of dis- 
order these modes may display a different degree of local- 
ization as, e.g., in a disordered photonic crystal (PhC). 32 
For a closed cavity these modes form a complete set and 
the field can be expanded in terms of the modes 



E = 



' N 



a n (t) exp (-iw„i)E„(r) 



(4) 



with E = {E x , E v , E z }. As far as a nonlinear po- 
larization is not present, the coefficients a n are time- 
independent. Conversely, in the general case, taking for 
Pjvl a standard third order expansion, one has for the 
non-linear interaction Hamiltonian 



H = -(J e V-P NL dV) 



-- — ^ !ft [Gjkim ajak^ a 7 . 

LOj+LO k =Ul+U>„ 



(5) 



where (...) is the time average over an optical cycle and 
the sum ranges over all distinct 4-ples for which the con- 
dition 



uij +uj k = L01+ uj v 



(6) 
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holds, with j,k,l,m = 1, . . . , N . The effective interaction 
occurring among mode-amplitudes reads: 



G 



jklrn 



-^/LOjLO k LL>iUJ m 
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(7) 
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with a,/3,7, S — x,y,z. This coefhcient represents the 
spatial overlap of the electromagnetic fields of the modes 
modulated by the non-linear susceptibility X ■ The dis- 
order is induced, e.g., by the random spatial distribution 
of the scatterers (as in random lasers) that leads to ran- 
domly distributed modes and, hence, to random suscep- 
tibilities and couplings among quadruple of modes. 

If the cavity is open, the mode set is no more complete, 
the modes whose profile decays exponentially out of the 
cavity are taken for the expansion all the others form 
the radiation modes. Under standard approach^— the 
coefficients in the expansion that weight the radiation 
modes can be expressed in terms of the disordered cavity 
one, and this results into linear terms in the Hamiltonian 
(open cavity regime). Thus, for an open cavity, Eq. ([5]) 
becomes 



^2 G fk a j a k + E G jkim a 3 a ka*ia* m 



j<k 



LUj -\-OJk — Ml +Cd T , 



(8) 



The Hamiltonian expressions, Eqs. ([5]), JSJ, can be 
also obtained starting from the corresponding Langevin 
dynamical equations, as detailed, e.g., in Ref. [Hj]: 

a n (t) = Y, G tn a 3+ E Gf^a^+rj^t) 

= + »,,«) P» 

where f] n (t) is a white noise, for which 

(r )j {t)rik(t')}=2T8 jk 8(t-t') (10) 



Here T is a "heat-bath" temperature, whose physical 
interpretation depends on the specific system. In the case 
of a random laser it represents the spontaneous emission 
and kgT = h/r, with r the amplifying level lifetimei 38 i 39 

Comparing Eq. © with the master equation for mode- 
locking lasers in ordered cavitie o 26 ' 33 

a n (t) = (g n - i n + iD n ) a n 

+{^-i5) ^2 a*a k ai +r] n (t) (11) 

01 j -\-OJ n = OJfc+LJ/ 

we can understand the physical role played by the param- 
eters of the probability distribution of the G's. Indeed, 
g n is the gain coefficient of the n-th mode in a round-trip 
through the cavity, £ n the loss term, D n the group veloc- 
ity of the wave packet, 7 the coefficient of the saturable 



absorber (responsible for passive mode-locking) and <5 the 
coefficient of the Kerr lens effect. Neglecting the latter we 
can see that a system with positive average of the Gj k i m 
corresponds to the presence of a saturable absorber. In 
the case of peaked probability distribution for the cou- 
plings P(G), i.e., weak disorder, the system will tend to 
display the same spectrum of many equally spaced modes 
typical of mode-locking lasers. In the present formalism 
this will be a ferromagnetic phase. One might, then, 
wonder what happens when the disorder is so strong to 
prevent the occurrence of this phase and, even, when the 
random coefficient corresponding of 7 is negative (i.e., 
when passive mode- locking is absent). We will discuss 
this issue in Sec. Ill Al 

In the "strong cavity limit", the linear coupling be- 

(2) 

tween modes is negligible and G m n is diagonal (i.e., one 
accounts only for the finite-life time of the modes) and 



' N 



iJ G 

=OJi~-\-LOi t 



(4) * * 



(12) 

Note that the modes in the disordered cavity may dis- 
play a different degree of localizations, as in the case of 
disordered PhC. Correspondingly, the distribution of the 
overlaps G spreads. Moreover, the constituents of the 
overlap integral are also very difficult to calculate from 
first principles. Indeed, to our knowledge, the only spe- 
cific form of the non-linear susceptibility has been com- 
puted by Lamt42i for a two-level system (without disor- 
der). Eventually, to estimate the coupling distribution 
from the experimental data is a very complicated inverse 
statistical problem, cf., e.g., Refs. 4l|f42| and references 
therein, and, so far, the reconstruction of the G's, for ex- 
ample, from measurements of random laser spectra has 
never been achieved. The interplay between susceptibil- 
ity and spatial distribution of modes leading to G's is, 
then, a very challenging problem that deserves a system- 
atic and sophisticated treatment that goes beyond the 
aim of the present work. 

In the following we will consider a mean-field approach 
in which all modes are connected among each other. We 
will, thus, approach the study of our model by means of 
Gaussian distributed G's with non-vanishing average, as 
detailed below. 

The leading regime considered in this work is, actually, 
driven by a quenched amplitude approximation, which is 
obtained by retaining the amplitudes A n = \a n \ (and cor- 
respondingly the energies of the modes) as slowly varying 
w.r.t. the phase <j> n — arg(a„), such that the resulting 
interaction Hamiltonian (retaining only those terms de- 
pending on the phases, and considering the strong or 
closed cavity regime, cf. Eq. ([T2j)). turns out to b o 27 i 38 

l,N 

W = — ^ Gi 1 i 2 i :3 i 4 Ai 1 Ai 2 Ai 3 Ai 4 x (13) 
cos(0 il + cj>i 3 - 4> l2 - faj 
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where the sum is limited to those terms that de- 
pend on the phases (i.e., we neglect terms whose indices 
are such that the argument of the cosine vanishes, e.g., 
ii = %2 and 13 = 24) and G is assumed real- valued. 

Actually, in the physical systems of our interest, it is 
not necessary that the resonant condition Eq. (|6|) for 
having four modes interact in the mode-locking regime is 
satisfied exactly. Indeed, it is enough that the mode com- 
bination tone Wi 1 lies inside an interval around Wi 2 +w^ 4 — 
u>i 3 corresponding to its linewidth4£ In the case, e.g., of 
the random laser, in which many modes oscillate in a rela- 
tive small bandwidth and are densely packed in frequency 
space so that the their linewidth overlap, this observa- 
tion supports the further mean-field-like approximation 
uji ~ wo, V i. In our model, therefore, the spectral dis- 
tribution of the angular frequencies will be considered as 
strongly peaked around wo and w^+w^ ~ Wi 2 +Wi 4 ~ 2wo 
so that the "selection rule" Eq. © is always satisfied. 

A suitable normalization and the introduction of an 
inverse temperature-like parameter j3 leads, eventually, 
from Eq. (T3D to 



PUjW\ = -P 



with 



E 



Ji cos(4> n + 4> i2 - (f> ia - cfo 4 ) (14) 
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Ji Ji 
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V 2 



(A 2 



(15) 

(16) 



© and ([16]). 

To the sake of qualitative comparison with the outcome 
of experiments the statistical mechanic inverse tempera- 
ture ft can be expressed in terms of the squared pumping 
rate as: 



V 2 
Jo 



(19) 



B. Finite temperature Bose-Einstein condensates 

A similar situation is found in the finite temperature 
Bose Einstein condensation with random potential. The 
zero temperature Gross-Pitaevskii equations 4 ^ reads as 



1T1 



<9$ 
~dt 



2m 



V 2 $ + F ext (r)$ + 5 |$| 2 $ (20) 



where V cx t(r) is an externally set disordered potential 
and g = Air£h 2 /m, with I being the s-wave scatter- 
ing length. An analogous model holds for reduced- 
dimensionality cases. 

The modes satisfy the time-independent linear 
Schroedinger equation 



— V 2 $„ + V oxt (r)$ n = E n <t> n 
Zm 



(21) 



Their interaction can be treated variationally by letting 



*(r,t) =^a„(*)$„(r)exp(-z^ 



(22) 



where Tbath is the heat-bath temperature, variance of 
the white noise 77, cf. Eq. (|10[) induced by spontaneous 
emission, and the squared volume factor guarantees ther- 
modynamic convergence (pH oc V ~ N). The average 
energy per mode is Eq = u>q(A 2 ). This is proportional 
to the so-called pumping rate V induced on the random 
laser by the pumping laser source. We will define it as: 



T> 2 - T ^ 

KrT\ 



2\2 



B -L bath 



(17) 



encoding the experimental evidence that decreasing the 
heat bath temperature^ 4 - or increasing the energy of the 
pumping light source 4 ^ has the same qualitative effect. 
The proportionality factor Jo in Eq. (|17[) is a material 
dependent parameter function of the angular frequency 
wo of the peak of the average spectrum, cf. Eq. ([8]), 



Jo = FWn 



d 3 r x (3) (^o;r) \E(r)\< 



(18) 



in which |i?(T")| = (E 2 ) ~ 1/V. Assuming that the non- 
linear susceptibility does not scale with the number of 
modes, the above integral scales as 1/V and Jq does not 
scale with the size of the system. The average of Ji, in- 
stead, scales as 1/V 3 , according to the definitions Eqs. 



A finite temperature model for BEC is the Stoof 
equation ) 47 ' 48 which is here written as 



l + h?fz K (r,t) 



2m 



V 2 $ + Fext(T-)$+ 5 |$| 2 $ 



(23) 



+ r){r,t) 



with fix = 1/fcsT (fcs is the Boltzmann constant) and 
where the finite temperature noise is such that 



(r,*(r>,t') V (r,t)) = ^ K (r, t)S(t - t')S^(r - r'). (24) 

Ti K {r,t) being the Keldish self-energy, which is imagi- 
nary valued (for its expression see Ref. [48j]) and nY, K cx 
~ l P~K ( see ^- e f- Ell)- Expanding over the complete set 
of the zero temperature equations, one obtains 



iha n (t) = —i^^aj n ajEje '» 



(25) 



+ ^2(G jk i n - iK jkln )ata j a k e-^ E ' +E "- E '- E ^ 
jki 

+Vn(t) 
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where rj n (t) — J d 3 r rj(r, t)<j> n (r, t), and the mode-overlap 
coefficients are defined as: 

Gjki m = g J $ 3 $ fc $i$ m d 3 r (26) 

and 

K jk im = J ^ K (r)^ k ^ m d 3 r. (27) 

Finally, the linear coupling coefficients come out to be 
a 3k = ^ f ^ K (r)^ k d 3 v (28) 

While retaining the synchronous terms (such that 
Ej + E k — Ei — E n = 0), the resulting equations are, 
hence, of the same form of those reported in Sec. II Al 
for the disordered electromagnetic cavity, being the en- 
ergy of the eigenstates in place of the angular frequency. 
Indeed, a strong coupling regime is attained when there 
is an enhanced region for the density of states. Con- 
versely, in other spectral regions, both the linear and the 
nonlinear coupling terms are averaged out by the rapidly 
oscillating exponential tails. 

Let us consider, for example, a periodic external po- 
tential with some degree of disorder. In this case, a Lif- 
shitz tail^ is present, that is, a region with energies in- 
side the forbidden gap corresponding to localized modes. 
This modes will all have approximately the same energy 
E = Eb where Eb is the band-edge energy, and will 
couple both among each other and with the delocalized 
Bloch modes at the band-edge. Correspondingly, the rel- 
evant equations for the strongly coupled modes are 

ifia n (t) = -i^2a jn ajE B (29) 

j 

+ ^2 (Gjkin - iKjUn) a*a,jak + r) n (t) 
jki 

The other modes (those far from the spectral gap) will 
be those mediating the thermal bath. The quenched 
amplitude approximation eventually leads to the phase- 
dependent Hamiltonian, Eq. (fT4")l . 

As discussed in the following section of the manuscript, 
even in the zero temperature limit a transition is ex- 
pected. This corresponds to the existence of a replica 
symmetry breaking transition in Bose Einstein conden- 
sates for finite and vanishing temperature, mediated by 
the degree of disorder and heuristically following the 
phase diagram reported in Fig. [T] below. 

C. Nonlinear optical propagation in disordered 
media and the zero temperature limit 

The nonlinear optical propagation of a light beam is 
described by the paraxial equation 



where A is the optical amplitude, k the wavenumber, n 
is the bulk refractive index and An is its perturbation 
due to disorder and optical nonlinearity (Kerr effect): 

An 

— = U(x,y)+n 2 \A\ 2 . (31) 

The nonlinear coefficient n 2 can be either positive (fo- 
cusing) or negative (defocusing) , while U(x,y) can be a 
perturbed (by disorder) periodical potential or a com- 
pletely disordered (speckle pattern) external potential. 
The resulting equation reads as 

r) A 1 n 
t-^A + U^A+^A-V. (32) 

This formally corresponds (with different meanings for 
the variable) to the zero-temperature two-dimensional 
limit of the Gross-Pitaevskii equations detailed above, 
cf. Ea. (|20]). 

In this well, the field can be expanded in terms 

of transversely localized (in two dimensions they are al- 
ways localized) electromagnetic modes, the energies be- 
ing replaced by their propagation wave-vectors. When 
there are bunch of modes such that their wave-vectors 
are similar, these will be strongly coupled and result into 
dynamical equations like Eqs. ©, (121)1) . This approach 
can be extended to three-dimensional propagation, en- 
compassing the dynamics of ultra-short pulses in random 
media as will be reported elsewhere. 

The replica symmetry breaking transitions investi- 
gated in the following will in general correspond to vary- 
ing coherence properties of the beam, eventually result- 
ing in unstable speckle patterns. The f3 — > oo limit of the 
statistical mechanical formulation of the problem has to 
be taken in this case (see, e.g., Ref. [5l[ for a simple 
case example in the framework of constraint satisfaction 
problems) . 



II. RANDOMNESS IN MODE-COUPLING 
COEFFICIENTS 

Let us consider our model Hamiltonian, Eq. (|13[) . in 
the mean-field fully connected approximation in which 
the non-vanishing components of the four index tensor 
Jii,i2,i3,ii — Ji are distributed as 

Ti = Jo/A 3 (33) 
(Ji-Ji) 2 = °*j/N 3 (34) 

The coefficient Jo was already introduced in the case of 
random lasers, cf. Eq. (fT8"|) . and A~ is the number of 
dynamic variables (mode phases) of the system, propor- 
tional to the volume V. The over bar denotes the average 
over the disorder. 

To quantify the amount of disorder, we introduce the 
"degree of disorder" parameter, i.e., a size independent 
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ratio between the standard deviation of the distribution 
of the coupling coefficients J\ and their mean: 

Rj = ?f (35) 

The limits Rj — > and Rj — > oo correspond, respec- 
tively, to the completely ordered and disordered case. 
The other relevant parameter for our investigation is the 
inverse temperature /3. For random lasers it is related to 
the normalized pumping threshold for ML, defined in our 
model as, cf. Eq. %HfyM 




where j3 = f3aj~± In general, /3 increases as the strength 
of nonlinearity increases or the amount of noise is re- 
duced. 



A. The ordered limit, saturable absorbers in 
random lasers, defocusing versus focusing 

With specific reference to the laser systems, as J 
grows the effect of disorder is moderated and for small 
enough Rj the model corresponds to the ordered case, 
previously detailed in Ref. }52j |. As also previously re- 
ported in Ref. [53j], a passive mode- locking (PML) tran- 
sition is predicted as a paramagnetic/ferromagnetic tran- 
sition occurs in /3. 

Indeed, in our units, when Rj — > 0, V = Vpml — 3.819 
(see Fig. [T]), in agreement with the ordered case^I — As 
explained below, the deviation from this value quantifies 
an increase of the standard ML threshold "Ppml due to 
disorder. The specific value for Vpml will depend on the 
class of lasers under consideration (e.g., a hber loop laser 
or a random laser with paint pigments) , but the trend of 
the passive ML threshold with the strength of disorder 
Rj in Fig. [T] has a universal character. The pumping 
rate V contains Jo : for a fixed disorder the threshold will 
depend on the nonlinear mode-coupling. 

A key point here is that the transition from continuous 
wave to passive mode-locking (PM — > FM) only occurs 
for a specific sign of the mean value of the coupling coef- 
ficient Jo, as shown in Fig. ^j. Comparing Eqs. © and 
(fTTj) one observes that this formally corresponds to the 
presence of a saturable absorber in the cavity (see also 
Ref. [26| and Sec. II A[) . In typical random lasers such 
a device is not present, and, hence, this ferromagnetic 
transition is not expected. 

On the other hand, the reported phase diagram, Fig. 
@ predicts that starting from a standard laser support- 
ing passive/mode- locking and increasing the disorder the 
second order transition acquires the character of a glass 
transition. A notable issue is that this phase-locking 
transition (normally ruled out for ordered lasers without 
a saturable transition), spontaneously occurs increasing 



/?, as an effect of the disorder and the resulting frustra- 
tion. 

With reference to nonlinear waves, the spontaneous 
phase- locking process is expected for a specific sign of the 
nonlinear susceptibility (corresponding to repulsive inter- 
actions for BEC and defocusing nonlinearities for optical 
spatial beams), for T = 0, amounting to Jq/cj > in 
Fig. (the threshold is at Jq/<jj = 4). For example, 
for a nonlinear optical beam propagating in a disordered 
medium, it is expected that above a certain degree of dis- 
order, there is a transition from a coherent regime to a 
"glassy coherent phase" , characterized by a strong vari- 
ation from shot to shot of the speckle pattern and, more 
in general, of the degree of spatial coherence. 



III. FUNDAMENTALS OF STATISTICAL 
MECHANICS OF DISORDERED SYSTEMS 

Hereby we report an extremely concise summary of 
ideas and techniques developed to deal with disordered 
systems. The aim is to let the non-expert reader find 
his/her way through the computation of the properties 
of our model that we present in Sec. IIVI and App. [X] 



1. Disorder and frustration: 
quenched disorder as technical tool. 

The main issue determining complex features, not 
present in ordered systems and involving collective pro- 
cesses that cannot be understood just looking at local 
properties, is frustration. This is usually a the conse- 
quence of disorder, not necessarily quenched disorder, 
though. Indeed, also in materials whose effective statis- 
tical mechanic representation is carried out through de- 
terministic potentials (as, e.g., for colloidal particles), a 
geometry-induced disorder can set up, determining frus- 
tration and a consequent multitude of degenerate stable 
and metastable states typical of glasses^— and spin- 
glasses^r— Quenched disorder, i.e., the explicit appear- 
ance of random coefficients in the Hamiltonian, allows an 
analytic computation, but the results are general and do 
not depend on the specific source of frustration. 



2. Statistical mechanics of a disordered system: 
the replica trick. 

In the presence of quenched disorder, one can compute 
the statistical mechanics of the system, averaging over 
the probability distribution of the disorder. In order to 
do this the so-called replica tricb^— can be adopted, or, 
else, the equivalent cavity methodj 68 ' 69 

The free energy of a single disordered system sample, 
denoted by J, is $j = — 1/ f3 log Zj. Correspondingly, the 
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physically relevant average free energy can be written as 



Complexity as a well-defined thermodynamic potential. 



log Zj 



where the overbar denotes the average over the distribu- 
tion of the J's. The latter coincides with the thermo- 
dynamic limit of any $j according to the self-averaging 
property required in order to have macroscopic repro- 
ducibility of experiments (the thermodynamics of a huge 
system does not depend on the local distribution of in- 
teraction couplings). 

To perform the average in Eq. (|37p is highly non triv- 
ial and one can proceed by considering n copies of the 
system, Eq. ([15]). 



(37) Besides numerous and hierarchically organized glob- 
ally stable states, glasses also display a large number of 
metastable states, that is, excited states of relatively long 
lifetime. In the mean-field theory such lifetime is, actu- 
ally, infinite in the thermodynamic limit because of the 
divergence of the free ener gy barriers with the size of 
the system, see, e.g., Ref. [73]. This means that, con- 
trarily to what happens in real glasses^ the number of 
metastable states at a given observation timescale does 
not change with time (after a given transient period) . Be- 
low a certain temperature (called dynamic or mode cou- 
pling temperature), the number Af of metastable states 
grows exponentially with the size ./V of the system (N 
being the number of modes in our cases). One can then 

(38) define an entropy-like function counting the metastable 
states as 



The average free energy per spin can, then, be computed 
in the replicated system, as 



(40) 



1 



lim — log Zj = 

JV-S-oo N 



Z 7 ) - 
lim lim — - 

N-yoan->0 N n 



1 



(39) 



where the average of the generic power of the partition 
function Zj is somehow computed for a finite integer n 
and, eventually, the analytic continuation to real n and 
the limit n — > arc performed. 



This is called configurational entropy in the framework of 
structural glasses, else complexity in spin-glass theory and 
its applications to constraint satisfaction and optimiza- 
tion problems. One can further look at the metastable 
states of equal free energy density /: Af(f) = expN'S(f) 
and at the free energy interval, above the equilibrium 
free energy / cq , in which the complexity is non zero: 

/e L/c q :/*]. 



3. Oddities of the replica formulation. 



Actually, to evaluate Z™ , one makes use of the saddle 
point approximation holding for large N (see Appendix 
IA1 for the specific case considered in this work) . That 
is, one practically inverts the limits N — > oo and n — > 
as expressed in Eq. (|3"9"1) . Yet, the method works. It 
took many years to rigorously overcome this oddity and 
a mathematical proof of the existence of the free energy 
can be found in Refs. I70ll7l1. 



IV. STATISTICAL MECHANICAL 
PROPERTIES 



Starting from the Hamiltonian, Eq. (|13p , replicated ac- 
cording to the prescription Eq. (1381) . and averaging over 
the disorder with the Gaussian probability expressed by 
Eqs. (|33| - p4|) . one obtains the following expression for 
the average of the n-th power of the partition function, 
cf. Appendix [A"l 



4- A probability distribution as an order parameter. 

The main novelty of the characterization of the spin- 
glass phase, historically first obtained by the replica 
method and subsequently confirmed by other methods, 
is that the order parameter is a whole probability distri- 
bution function describing how different thermodynamic 
states are correlated. The degree of the correlation be- 
tween two states is called overlap. In mean-field the- 
ory different states exist that can be more or less corre- 
lated according to their distance on a tree-like hierarchi- 
cal space called ultrametric^ 



Z n j = JvQVA e~ N n G[Q < A] 
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(41) 



n G [Q, A] = n A [Q, A] + log Z^ [A] 
n A [Q, A] = £(l + \r a f) - PI* £ \fh a \ 



a=l 
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A j 'ft J. j ft 

^2{ltb+ l r a&| 4 ) - [<labKb + 5ft {r ab p, ab 



a<b 



a<b 



+ m a v a \ 



(42) 



a=l 
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Z+[A] = f n# e-^*[«>> A ] (43) 
J 0=1 

l,n 

-/8W aff [{^};A] = 5^K[e^--^A o6 + e t W« + *»^ a6 ] 
+ ^5i[e 2 ^ a + e^ a ] (44) 



Eq. (061): 



a=l 



2?Q = fj N 2 dq ab dr ab x J| N 2 df a drh a 

a<b a— 1 

where Q = {5, r, f , to} and A = {A, fi, ft, v\. The overlap 
matrices q ab , X a b are real- valued, whereas the others have 
complex elements. 

The integral Eq. (|4Tj) is evaluated by means of the sad- 
dle point approximation (valid for large N). The above 
expressions need a form of the matrices q abl r abl X ab and 
fi ab to be completed. Contrarily to what might seem 
reasonable, the form providing the thermodynamically 
stable solution is not the one in which all replicas are 
equivalent, i.e., all elements in the matrices q ab , r ab , X ab , 
li ab are equal. One must, thus, resort to a spontaneous 
Replica Symmetry Breaking. In Appendix [X] we report 
the computation of thermodynamics both in the Replica 
Symmetric (RS) approximation and in the "one step" 
Replica Symmetry Breaking Ansatz (1RSB), i.e., the ex- 
act solution for the system under probe. In the following 
we, thus, analyze the properties of the latter solution. 

Spin-glass systems described by more-than-two-body 
interactions, cf. Eq. (|13p . are known to have low 
temperature phases that are stable under the 1RSB 
Ansatz . 75 i 7692 Under this Ansatz, taking the n —> limit, 
the free energy functional /?$ reads, cf. Appendix lAl 

/3$(m;Qg),Ag)) = G(m; Q£>,A&>) = (45) 



1 ~ ,4 
•— — m 
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-m( go 4 + |r | 4 ) 



l-(l-m) (^ + |rr| 4 ) 
1 — m 



/ \ — \ — — ~ 

+ y (A <?o + ^o?"o) - V-aXd ~ vm 



(\\q\ + nm) 

Ai 
2 



J 2?[0]log y jf # exp£(0;O,l) 

where = {x , C(f , Co }, 1 = {^i,Cf ,Cf}, 2?[o] is the 
product of three Normal distributions and 

) = «{e* [Ci VAA - |A M | + Co v/Ao-l^ol + 

t)} <«» 



£(</>; 0,1) 



with AA = Ai — Ao, A/i = — /iq. For later convenience 
we define the following averages over the action e c , cf. 



Jq 2 " dcf> cos <fi e *~ 
J„ 2lr d</> sin 



cc 



sc 



(sin < 



— Jo 



(47) 
(48) 



The values of the order parameters Ao,i, (J-o.i, fJ-d and v 
are yielded by 



Ao.i = — (<?o,i) 



M0.1 



/3 2 



ko,i| 2 ^o,i (49) 



s 2 



\r\ 2 f 



/3Rj 1 - ,2 ~ 
— - — to m 



(50) 



The parameter m (without tilde!), whose meaning will be 
discussed below, takes values in the interval [0, 1]. The 
remaining parameters are obtained by solving the self- 
consistency equations: 





= «4>m)o + 


«4 


m)o 




(51) 


Qo 


= ((C£) 2 „) + 


((sc 


m)o 




(52) 


1*1 


= «4)m)o - 


«4 


m)o + 


2«((c£S£) m )o 


(53) 


ro 


= ((cc)l)o 


((sc 


m)o "+ 


M( c £)m)o((sc) 


m)o (54) 


f - 


= (((e^cU) 


0; 


m = 


({{e z *)c)m)o 


(55) 



where the averages are defined as 
<(..•)>« = 

((...))o = / 2>[0](... 



/©[!](..•) 




m 






m 



(56) 
(57) 



These equation are solved numerically by an iterative 
method. The overlap parameters qo,i are real- valued, 
whereas ro,i,f and to are complex. "One step" parame- 
ters Xo,i (A — q, r) enter with a probability distribution 
that can be parametrized by the so-called replica sym- 
metry breaking parameter to, such that 

P(X) = to 5(X - A ) + (1 - m)5(X - Ai). (58) 

The resulting independent parameters (there are ten of 
them) that can be evaluated by solving Eqs. (|5ip - ([55|) 
must be combined with a further equation for the pa- 
rameter to. This is strictly linked to the expression for 
the complexity function of the system. 



COMPLEXITY 



In the order parameter Eqs. (j4"Tj)) - (f55|) to is left unde- 
termined. An additional condition is needed to fix the 
value for this parameter. The first possibility is treating 
to as a standard order parameter: in this case the thermo- 
dynamic state corresponds to extremizing the replicated 
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free energy (thus, maximizing it^), i.e., implementing 
the self-consistency equation 



d$(m;Q sp , A sp ) 
dm 







(59) 



The highest temperature at which a solution exists with 
to < 1 furnishes a transition temperature between para- 
magnet and glassy phase: the Kauzmann or static tem- 
perature (T s ). This is an equilibrium thermodynamic 
phase transition. 

This approach, however, does not reflect the known 
physical circumstance that a glassy system exhibits ex- 
cited metastable states also at temperature T above T s j2£ 
where the equilibrium phase is paramagnetic. Vitrifica- 
tion, indeed, is due to the presence of a not vanishing 
complexity at a temperature above T s (and below some 
Td > T s ), i.e., to the presence of a number of energetically 
equivalent states with free energy / > / eq (T). Since, 
however, energy barriers tend to infinity in the thermo- 
dynamic limit in the mean-field approximation, the sys- 
tem dynamics is forever trapped in one of these states 
for T < Td- The temperature Td is, thus, called dynamic 
transition temperature. 

Across this transition the complexity E, defined in 
Eq. pTJ)) . starts being different from zero. Exactly at 
T = Td the complexity as a function of free energy, £(/), 
has a delta-shaped non-zero peak at the free energy /i 
which corresponds to a maximum of S(m) for a value of 
in = m(fi) = 1. In our 1RSB formalism: 



9S(m;Q S p, A sp ) 
dm 







(60) 



As T decreases (T s < T < Td), the complexity is not van- 
ishing for an increasing range of free energies f* > f > fi 
that corresponds to a range for to: m* < m < 1. The 
complexity shows a maximum for to < 1 at to* (/ = /*) 
solution of dYijdm = 0, while it is at its minimum value 
for to — 1 and f = fi- We stress that this is not a solu- 
tion to Eq. ([59"]) . 

Lowering the temperature, at T = T s the minimum 
value of complexity - corresponding to to = 1 - vanishes, 
i.e., it is a solution to Eq. (|59|) . and /i = / eq corresponds 
to the free energy density of the global glassy minima of 
the free energy landscape: as mentioned above, we are 
in presence of a thermodynamic phase transition and the 
thermodynamic stable phase is a glass. 

The physically significant value for to is to*, corre- 
sponding to the maximum of S. It denotes the value of 
free energy /* where the number of states is maximum 
and exponentially higher than the number of states at 
any / < /*, and, hence, the most probable (among those 
of the metastable states). At the thermodynamic tran- 
sition point from the paramagnetic state to the glassy 
(T = T s ) it holds f PM = fi = /eq = *- 95 Below T s 
fi < f cq (hence, £(/*) < S(/ eq ) = 0) and the physically 
relevant £(/) has a support [/ e q,/*]- 

In the following we will analyze the whole complexity 



vs. free energy curve £(/) at given /3, Jo and the behav- 
ior of the minimal positive complexity £(T) (and S(7 5 )) 
between T s and Td- 



A. Computing the complexity functional 



In Eq. (|40|) one needs to know the number of 
metastable states, that are the local minima of the 
free energy landscape. Would we know the landscape, 
though, we would have solved the problem already. If 
self-consistency equations for local order parameters are 
known, a possible analytic approach to get information 
on the complex landscape is to guess a trial free energy 
functional whose stationary equations lead back to the 
self-consistency equations. This is what Thouless, An- 
derson and Palmer (TAP) proposed in the framework of 
spin-glasses starting from the self-consistency equations 
for local magnetizations. 77 Starting from TAP functional 
and TAP equations and considering solutions to the TAP 
eqs. as states (with some assumptions to be a posteriori 
satisfied) one can build the functional £ from Eq. (|40p . 
cf., e.g., Refs. (zHH- 

A comparative study to the TAP-derived complexity 
functional and the replicated free energy, computed in a 
general scheme that includes the Parisi Ansatz, 86 allows 
to show that the Legendre Transform of <E> with respect 
to the single state free energy coincides with Eq. P0"j) . 
According to this approach, in our model the complexity 
can, thus, be explicitly computed as the Legendre trans- 
form of Eq. (j45)) : 



(61) 



E(to; Q sp , A sp ) 

= min [— /3to<I>(to) + /3to/] 

TCI 

a 2&* 

= /3m — 
am 

= ^2 m2(|5i| 4 + |ri| 4_ ko |4_ |ro |4 ) 



P[0]log J V[l] dtp exp £((/>; 0,1) 

P[0](log / dej) cxp£(0;O,l)) m 
Jo 



where the single state free energy 

<9(m$) 



/ 



dm 



(62) 



is conjugated to to. Since the above expression is propor- 
tional to d^/dm, equating S = provides the missing 
equation to determine the order parameters values. 



VI. PHASE DIAGRAM AND COMPLEXITY 

By varying the normalized pumping rate V and the 
degree of disorder Rj, we find three different phases, as 
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FIG. 1: Phase diagram in the V,Rj plane. Three phases are 
present: PM (low V), FM (high V /weak disorder) and SG 
(high "P/strong disorder). The full lines are thermodynamic 
transitions, the dashed line represents the dynamic PM/SG 
transition. 
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FIG. 2: Phase diagram in the plane Jo,T in aj units. Also 
negative Jo are considered. Three phases are found: PM 
(high T, low Jo), FM (low T/large J ) and SG (low T/low or 
negative Jo). The full lines are thermodynamic transitions: 
random first order between PM and SG and standard first 
order between PM and FM and between SG and FM. The 
dashed line represents the dynamic PM/SG transition. 



shown in Fig. [T]in the ("P, Rj) plane and in Figs. [3] in 
the (T, J ) plane. 

Paramagnetic phase — For low V the only phase 
present is completely disordered: all order parameters are 
zero and we have a "paramagnet" (PM); for the random 
laser case this phase is expected to correspond to a noisy 
continuous wave emission, and all the mode-phases are 
uncorrelated. Actually, this phase exists for any degree 
of disorder and pumping, yet it becomes thermodynami- 
cally sub-dominant as V (or /3) increases and, depending 
on the degree of disorder, the spin-glass or the ferromag- 
netic phases take over. 

Glassy phase — For large disorder, as V / (3 grows, a dis- 



FIG. 3: (Color online) Detail of the Jo,T phase diagram 
around the tricritical point. Full lines are thermodynamic 
transitions. Also the transition between the SG (1RSB) 
and the approximated RS solution for the FM phase is dis- 
played (double-dotted line) showing no appreciable difference 
with the exact one. The dashed line represents the dynamic 
PM/SG transition. The dotted bold line represents the FM 
spinodal lines both inside the PM and the SG phases. The 
spinodal of the RS FM phase is shown as well (smaller dots). 



continuous transition occurs from the PM to a spin-glass 
(SG) phase in which the phases (f> are frozen but do not 
display any ordered pattern in space. First, along the line 
Vd = y/0d/Rj, in Fig. m or at T/aj = l/p d = 0.15447 
in Figs. [2] [3] (dashed lines) a dynamic transition oc- 
curs. Indeed, the lifetime of metastable states is infinite 
in the mean-field model and the dynamics gets stuck in 
the highest lying excited states. The thermodynamic 
state is, however, still PM. Fig. [3] displays a detail of 
the tricritical region where. There, besides thermody- 
namic transition lines, we also plot as dotted curves the 
lines at which the ferromagnetic phase first appears as 
metastable, i.e., the spinodal lines. 

In Fig. H] we plot the complexity of the metastable 
glassy states of lowest free energy between the dynamic 
and the static transition. In the left panel of Fig. [5] 
is displayed for three different values of Rj; the thresh- 
old pumping for non-zero minimal complexity grows as 
the degree of disorder Rj decreases, as well as the cor- 
responding V range. In the right panel £(T) is plotted 
and it is independent of Rj. In Figs. [6] and [7] we display 
two instances of the whole complexity curve both vs. / 
and m at T — T s and at a higher temperature T < T d . 

Across the full line V s {Rj) =js/K/Rj, in Fig. fflor, 
alternatively, across T/aj = l/fl s = 0.14099 in Fig. H 
a true thermodynamic phase transition from the contin- 
uous wave (paramagnetic) phase to the "glassy coherent 
light" (spin-glass) phase occurs. The order parameter 
qi (the Edwards- Anderson parameter qea— ), discontin- 
uously jumps at the transition from zero q\ > qo = 0, 
while fh = r = rj = r d = (see Fig. [51 bottom panel). 
The SG phase exists for any value of Rj and /3 > j3 s . 
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FIG. 4: Complexity E('P) of the lowest lying glassy states 
in free energy between the values of the pumping rate corre- 
sponding to the dynamic and static transition from the PM 
to the SG phase at Rj = 0.5. 
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FIG. 5: Left: Complexity of the lowest lying glassy states in 
free energy E('P) between dynamic and static transition from 
the PM to the SG phase along Rj = 0.3, 0.4 and 0.5 lines. 
The qualitative behavior is identical for any Rj > 0.3. Right: 
£ vs. the effective temperature T in aj units. 
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FIG. 6: £(/) (left) and E(m) (right) in the glassy phase at 
the static transition effective temperature, T = 0.14099. This 
is the picture holding for any Rj > 0.3. 



FIG. 7: £(/) (left) and £(m) (right) in the glassy phase at the 
static transition effective temperature, T = 0.15 < Td- The 
lowest state free energy of metastable glassy states is denoted 
by /i (i.e., corresponding to m = 1 in the right panel, see 
text). The free energy of maximum complexity is denoted by 
/*, correspondingly m* in the right hand side plot. 



In the stable SG phase, metastable states (with infinite 
lifetime) continue to exist so that the thermodynamic 
state is actually unreachable along a standard dynamics 
starting from random initial condition. In Fig. [3] we plot 
the typical behavior of the complexity versus the single 
state free energy at T/aj = 0.1, qualitatively identical 
to the left panel of Fig. ©displaying £(/) at T = T s . 

Ferromagnetic phase — For weak disorder a ran- 
dom ferromagnetic (FM) phase turns out to dominate 
over both the SG and the PM phases. The transition 
PM/FM line is the standard passive ML threshold (see 
e.g. [HUll) and it turns out to be first order in the 
Ehrenfest (i.e., thermodynamic) sens e 74 ' 89 . From Fig. 
Q] we see that it takes place at growing pumping rates 
V for increasing Rj until it reaches the tricritical point 
with the SG phase. In the (T, Jo) plane it occurs at large 
- positive - Jo, cf. Figj2] 

To precisely describe the FM phase in the 1RSB Ansatz 
we have to solve eleven coupled integral equations [Eqs. 

and Eq. J55|) (E(m;Q sp ),A sp ) = 0), cf. Eq. 
([61]) ]. In evaluating their solutions we have to consider 
that, in the region where the FM phase is thermodynam- 
ically dominant, both the PM and the SG solutions also 
satisfy the same set of equations. Besides, unfortunately, 
the basin of attraction of the latter two phases - in terms 
of initial conditions - is much broader than the FM one. 
Starting the iterative resolution from random initial con- 
ditions, determining the FM transition and spinodal lines 
becomes, thus, numerically demanding. 

An approximation can be obtained by considering 
the Replica Symmetric (RS) solution for the FM phase 
(FM rs ). This reduces the number of independent param- 
eters to seven (qi = qo, r l ' = r Q : , r d ' and rh 11 ' 1 ). The 
corresponding transition line is shown as a 
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FIG. 8: Discontinuity of the order parameters at the transition points for three values of Rj. Top Left panel: jump in qo,i, 
at the PM/FM transition in V for small disorder, Rj ~ 0.1; top right: discontinuities in ro,i, f and m at the same transition. 
For such small Rj the replica symmetry breaking is practically invisible: q\ ~ go, r\ ~ ro [to the precision of our computation, 
O(10~ a )]. Mid left panel (across tricritical region in Fig. [TJ qo.i vs. V at Rj ~ 0.26 where, increasing the pumping rate, first 
a PM/FM transition occurs followed by a FM/SG one. Mid right panel: ro,i, f and m vs. V for the same interval. First order 
transition point are signaled by vertical lines. Left bottom panel: go,i vs. V for large disorder, Rj = 0.4 across the PM/SG 
random first order transition. Right bottom: ro,ii r and rh are always zero in the SG and in the PM phase. 
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dashed-dotted line in Fig. [3J where, around the tran- 
sition, we observe no practical difference with the exact 
SG/FM, even though the replica symmetry is clearly bro- 
ken. 

In Fig. [8] we show the discontinuous behavior of the 
order parameters across various transitions. As disorder 
is small (top panel) one can observe that the RSB of 
the solution representing the passive mode-locking phase 
vanishes, at least for what concerns the limit of precision 
of our computation. As the degree of disorder takes val- 
ues around the tricritical point the RSB is clearly visible 
(mid panel), both in the FM and in the SG phases. For 
increasing disorder the FM is absent (Rj > 0.263) and at 
high pumping/low temperature only the glassy random 
laser phase remains. 

We must necessarily implement the 1RSB Ansatz, 
though, to determine the not-vanishing extensive com- 
plexity which signals the presence of a large quantity of 
excited states with respect to ground states and study 
its behavior in T and Rj. This, as anticipated, also im- 
plies the occurrence of a dynamic transition besides the 
thermodynamic one. In the phase diagrams, Figs. [TJ [21 
[31 this takes place between PM and SG, where the state 
structure always displays a non-trivial S, for any (3 > fid- 
Whether an exclusively dynamic transition can occur as 
a precursor to the FM phase, as well, could not be di- 




0.02 0.04 0.06 Ht, 0.0005 0.001 f - f eq 



FIG. 9: Complexity vs. free energy curve is plotted in the SG 
phase (left) at T/aj — 0.1 and in the FM phase (right). 



rectly established in the present work. Indeed, the region 
of expected dynamic transition lies beyond the spinodal 
FM line, already very difficult to obtain numerically be- 
cause of the competition with the SG and PM solutions. 
However, the existence of a metastable FM phase (cf. 
spinodal line in Fig. [3j with an extensive complexity, cf. 
e.g., Figs. [51 and [TU1 might well correspond to an arrest 
of the dynamic relaxation towards equilibrium of the sys- 
tem. 

In the right inset of Fig. [S]we show, e.g., £(/) in the 
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FIG. 10: Complexity curves of the FM phase at Rj = 3.54 
at temperatures between T = 0.082<jj (right most) and 
T — 0.139<t j (left most). Both the magnitude of the maximal 
complexity and the free energy interval in which £(/) > 
decrease. Notice that the equilibrium free energy decreases 
as temperature increases. 
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(left inset of Fig. [9]) that is sensitively larger and does 
not depend on the Rj: the maximum complexity drops 
of about two orders of magnitude at the SG/FM transi- 
tion, thus unveiling a corresponding high to low complex- 
ity transition. 

In Fig. QTJ at a relatively low temperature T — 0.0785 
we show the behavior of the 1RSB (equilibrium) free en- 
ergy and order parameters across this SG (high complex- 
ity) /FM (low complexity) transition. The transition is 
first order in Rj. 



VII. CONCLUSION 



We have reported on an extensive theoretical treat- 
ment of the thermodynamic and dynamic phases of non- 
linear waves in a random systems. The approach allows 
to treat nonlinearity and an arbitrary degree of disorder 
on the same ground, and predict the existence of complex 
coherent phases detailed in a specific phase-diagram. The 
whole theoretical treatment is limited to the quenched- 
amplitude approximation, which allows to catch the ba- 
sic phenomenology and to demonstrate the existence of 
phases with a not-vanishing complexity in a variety of 
physical systems, and specifically random lasers, finite 
temperature BEC and nonlinear optics. This approxi- 
mation will be removed in future works, and novel exotic 
phases of light in nonlinear random system will be de- 
tailed. 

Our theoretical work shows that the interplay of non- 
linearity and disorder leads to the prediction of substan- 
tially innovative physical effects, which bridge the gap 
between fundamental mathematical models of statistical 
mechanics and nonlinear waves. This allows to identify 
frustration and complexity as the leading mechanisms for 
a coherent wave regime in nonlinear disordered systems. 
Natural extension of this work will be considering the 
quantum counterpart of the predicted transitions, and 
the analysis of out of equilibrium nonlinear waves dy- 
namics. 
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FIG. 11: Top panel: free energy of the FM and SG phases 
vs. Rj at T — 0.0785. As the degree of disorder increases 
the system undergoes a first order phase transition from a 
ferromagnetic phase to a spin-glass. Mid panel: order param- 
eters r*i,ro, r and the magnetization m are shown vs. R.j. 
Beyond the transition point their values drop to zero in the 
SG phase. Bottom panel: q order parameters for the FM and 
the SG phase. 

FM phase at (Rj,V) = (0.28, 5.92). This has to be com- 
pared with the SG complexity at the same temperature 
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Appendix A: Replica computation of the 
thermodynamic properties 



The replicated partition function of the system de- 
scribed by the Hamiltonian H[{(f>j}}, cf. Eq. (p~3|) . reads 



ryri 

A/ — 



n N 

nn< 

a=l J = l 



exp 



-0 



0=1 



(Al) 



In order to compute the free energy of the system using 
the replica trick, cf. Eq. ([551) . Eq. (jAll) has to be av- 
eraged over the probability distribution of i.i.d. random 
bonds: 



P(J) 



I N 3 



exp 



-JV 



3 (J- J a /N^ 



2aj 



(A2) 



Eq. (|4"T1) . then reads 



X>0expjiV 
16 



^t( 1+ |*.««)f) 

a— 1 x 7 
E(l^({0})| 4 + |i2a6(W)r 



0=1 



(A3) 



with 



JV 



(A4) 



a=lj = l 



and 



JV JV 
Q ab {{<P}) = -J2e<f?-^; M a {{<p}) = ^ 5>*? 



(A5) 



1 JV 1 JV 



(A6) 



where we used the Euler's formula to represent the co- 
sine and introduced abbreviations for the quantities Eqs. 
(|A5|) - (|A6|) . We notice that the matrix Q ab is Hermitian. 
A further step is to introduce extra parameters - that will 
eventually result as the order parameters identifying the 
various phases of the system - by means of the following 



identities: 



1 = II N I [N {Qab({(f>}) - qab)} (A7) 

a<b •* 

l = tl N dTabS i N i R M}) ~ r ab )} (A8) 

a<b J 

n „ 

1 = 11^ dr a 6[N(R a ({<t>})-r a ^] (A9) 
0=1 

n . 

1 = II N / dThJ \- N ( M «(M) - "*<»)] (A10) 
0=1 J 

5{N(Q ab ({<p})-q ab )} = 



5[N(R ab ({<j)})~r ab )] 
N (R a ({<t>}) - r a 



dXab JH[K b N(Q ab (U})-Qg b )} 

(All) 
(A12) 

d ^ jt[n„N(R a (u})-fg)} 



2tt 



(A13) 



d[N(M a ({<f>})-m a )} = I 



dv a „SRfa 



SR[C al> JV(M ({fli})-ma)] 



The two-index auxiliary variables q ab , r ab , defined for dis- 
tinct couples of replicas a and b, with a < b, will be 
considered in the following as elements of symmetric ma- 
trices [cf. Eqs. (|A26j) and (|A52ft ], i.e., g a fc = q ba and 
^ab = r ba . In particular, since Q a b{4>) = Qba(<f>), Eq. 
(jA7|) implies that q ab = q ab are real valued. 

Denoting for shortness the sets of parameters by the 
"vectors" Q = {q, r, f, m} and A = {A, fi 7 jl, u}, this leads 
to 



Z n j = JvQVA e~ N n G[Q < A] 

n G [Q, A] = n A [Q, A] + log [A] 



(A14) 



n A [Q, A] 



2^2 " 



32 



if„i 4 



a=l 



o=l 



16 



l,n 



'' X ' ] ~\b\ 4 + \r a b\ 4 ) ~ ^ ^[gabAqj, + r ab \x ah \ 



a<b 



a<b 



(A15) 



-}^\fgfj, a +fh a V a ] 

a=l 

/n 
J] d^e-W^W (A16) 
0=1 



.^5R[e 2 ^A a + e^^] (A17) 
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with 



and this trivially leads to the identification 



l,n 



VQ = Y[ N 2 dq ab dr ab x [| N 2 df a dm a 



a<b 
l,n 



jyj^ _ rr dX a b dfi a b ^ y-r djl a dv a 

= 11 OlT 9tT X 11 



a<b 



27T 2tt 



The average replicated partition function integral can be 
estimated by the saddle point method for large N, i.e., 
by approximating 



X>[X] e- 

Y Y =0 V j = 1, . . . , # parameters 



ax, 



Denoting by (. . .) e g the average over the measure e j 
the saddle point equations are: 



4 



9ab| <?ab (A18) 

r a& | 2 r a , (A19) 



V a, b = 1, . . . , n; a < 6 



(e 2 ^) eff ; 



Aa = ^|r„| 2 f a (A20) 
o 

= \rh a \ 2 m a (A21) 



V a = 1, . . . , n 



The diagonal values of the overlap matrices are set to 
zero. Eventually, according to Eq. (159"1) . one has 



I e -nJVG[Q ap ,A sp ] _ 1 

/3$ = -— lim = lim G [Q sp , A sp ] 

(A22) 

The parameters with a single replica index turn out not 
to depend on the specific replica. Indeed, Eqs. (|A20[1 - 
(|A2ip might in principle be obtained by perturbing the 
original Hamiltonian with a small field coupled to a local 
function of the planar, XY, spins S — e 1 ^, independently 
from the possible introduction of replicas. 



r<2 



If the perturbation is — fej^ - =1 S 2 we obtain 

s ] -> (8) ( as N -> 00 ) 



1 aiogZj 
A 9fe 



A ^— ' J 



(A23) 



1 n 

r a = lim - Y] r a 



(A24) 



a=l 



Similarly, perturbing Eq. (fT3")) with — ^X)j=i we S e * 

1 SlogZj 1 A — 
m = — tt, — = — y Si (A25) 
N dh N i— 1 y ' 

19^ If. 

m = — lim — — = lim — > m a 

n->o nA oh n->o n ^— ' 

a— 1 

Though no external ad hoc perturbation can be ap- 
plied to the Hamiltonian Eq. (|A17I) to reproduce two in- 
dices quantities, the same symmetry should apply, since 
all replicas of the original problem were introduced in 
the same way: the system is symmetric under replica ex- 
change. This is called the replica symmetric (RS) Ansatz. 

Qab = q Va ^ b ; r ab = r \/a^b (A26) 



1. Replica Symmetric Ansatz 
In this Ansatz Eqs. (|A16IA17|) become 



« n 

zf = j n 



/3H c g = — | n 



-/3« rff [{<M] 



-5ft 



E< 

a=l 
2 



,l<f>a 



(A27) 
(A28) 



The second term in the rhs can be rewritten as 



= 5ft 



(A29) 



:■)? 



\ a=l a=l / 



2 



E< 

o=l 



that is valid for any replica, and, therefore independent 
from any replica index: f<j = f a , Va = 1, . . . , n. In the 
replica formalism, the same quantity can equivalently be 
written, as 



lim 



1 dZ n , 



n ^ 

3 - = lim - } r a 



o nN dk n->o n 



The squared terms in the exponent of the integrand can 
be linearized by using 



>| 2 /2 _ / dC R dC J „-\C\ 2 /2„dXCw) 



2?r 



2tt 



(A30) 
(A31) 
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thus yielding 

zf = J v P {x)Vp{Q R )Vp{?) J 

C(<p;x,() = 3? 



(A32) 



\n\+Xy/2fL + v) (A33) 



2? (t 



(*-§) 



Vp(w) 



dw 



The replicated free energy eventually reads: 
/3$ 



■ _|?| 4 -|r| 4 + |if| 



32 



(A34) 



(A35) 



■ ^|m| 4 + ^(1 - - ^ - 2£f - 2Dm] 



Vp{x)Vp(C, R )Vp{C, 1 )\og / # e £W;;r ' C) 



Deriving w.r.t. to Q's parameter we obtain the specifi- 
cation of Eqs. (|A18HA2lj) for the replica overlap param- 
eters q, r and for f and m 



(A36) 
(A37) 



A = 






A* = 


4 |r 


2 

r 


A = 


8 |r 


2 - 

r 


v = 


)8Jb - - , 

TO 

2 1 1 


2 - 
TO 



Taking the derivative of /3<& in Eq. (|A35I) w.r.t. Ji and 
v we obtain 



r d = ((e 2 ^) £ ),, c 
m = ((e l0 >£),, C 



where we define 



r 2ir 



<,C(<t>;x,C) 



(A38) 
(A39) 



(A40) 



Deriving G w.r.t. A and /i and equating to zero we 
obtain 

Q R = (cl+sl) xX (A41) 

r = (4- 4 + 2* c £S£ ) x C (A42) 

C£ = (cos0) £ S£ = (sin</))£ 

after having integrated by part in the Gaussian measures. 
To help the non-expert reader to easily derive the self- 
consistency equations we exemplify the calculation of Eq. 
(IAT81) . 



dG 



dX R 



= = 1- q ] 



(A43) 



The latter term can be simplified by integrating by part 

dF(y) 



vp(y) y F[v) = / Vp(y) 



dy 



(A44) 



with y = C R , C 1 in Eq. (|A44j) . yielding 

((V+64c=V^~W x (A45) 



(cos 2 0-c £ +sm 2 (/>~s 2 c ) x 



C 



The self-consistency equation can thus be rewritten as, 
cf. Eq. ((ST]), 



l-^ = l-(4 + 4) 



We recall that since q a b is real, and so is X a b, cf. Eq. 
(|A18p . in the RS Ansatz the equations q 1 = X 1 = 0. 

Before deriving Eq. (|A42I) , we rewrite the part of Eq. 
(IA33[) involving the integrating variable x as: 



(A46) 



xv^cos^l + g+sm^l-^j 

In determining the above expression one can use, e.g., 
the trigonometric law of tangents to yield 



2 /i + 



(A47) 



and the relationships between trigonometric and inverse 
trigonometric functions: 

a 

sin[arctan(6')] 



cos[arctan(0)] 



1 



Using Eq. (|A47j) . together with Eqs. (|A38|) and (|A41(1 . 

we have: 



<9G 



9/i 



; ,-0 = ^-^ + ^(l-g) 



a; c £ W 1 



(A48) 



Integrating by part with Eq. (|A44|) . y = x, we find 



A Cc f + V\ +Sc f-¥\, 



(A49) 



(cos 20) c - c 2 c + s 2 c 

R 

+ ^ (i-4-4) 
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and eventually one obtains the real part of Eq. (|A42I) . 
The imaginary part of the self-consistency equation for r 
is analogously determined from Jp- = 0. 

Above a given critical temperature (depending on Jo) 
the solution to Eqs. (|A38|A39|A41|A42|) is paramagnetic, 
i.e., q = r = rh = f = and the free energy is 



PM 



32 



l0g27T 



(A50) 



Below T c (Jq), depending on the value of Jo the solution 
can either be ferromagnetic m / or spin-glass to = 
0. The latter solutions are, however, not stable against 
fluctuations in the space of replica overlaps? 6 and, thus, 
we have to try an Ansatz different from Eq. (|A26[) to 
provide a self-consistent thermodynamics. 



2. One step of Replica Symmetry Breaking 

In order to obtain a thcrmodynamically consistent re- 
sult the symmetry cannot be conserved. We are in pres- 
ence of a spontaneous Replica Symmetry Breaking. The 
way to break the symmetry must be a-priori hypothe- 
sized, since there has been found, so far, no way to deduce 
it. The correct way to express the elements of the over- 
lap matrices is called Parisi Ansat a 66 i 67 and, depending 
on the kind of system, can consist of one or more RSB's. 
According to what happens in other spin models with p- 
body quenched random interactions (p being larger than 
2), the right Ansatz for the matrices of our model is the 
one-step RSB, that is, we have a n x n matrix divided in 
square blocks of to x to elements 

q a b = qi, r ab =n if/(— )=/( — ) (A51) 

Vto/ \to/ 

Qab = <7o; r ab = r if / (—) ^ I ( — ] (A52) 

Vto/ \ m / 

For instance, for n = 6 and to = 3. 



/ qi qi q q q \ 

qi q\ qa qa qa 

qi Qi q Q q q Q 

qa qa qa qi qi 

. qa qa qa qi <?i 

V 9o <7o qa 9i ?i / 



?(a/3) 



The one replica index observables are instead still RS, 
as exemplified in Eqs. (IA23IIA25]) . Now, let us write the 
"vectorial" replica index 

a -» a = (oi, a 2 ) 



n n/m rn 



a=l 



ai — 1 a2 — 1 



Take a 1RSB matrix K ab and two replicated observables 
g a and h a . The following expressions hold for the sum of 



a generic product 

2 K ab g a h b = K ab g a h b 



(A53) 



a<b 



n/m l,m l.n/m l,m 

= Kl ^ ^ ga 1 a 2 ha 1 b 2 + Ka ^ 9a 1 a 2 hb 1 b 
ai=l a 2 ^-b 2 aiT^bx a 2 b 2 

n/m l,m n/m m 

K\ ^ ' ^ ' 9a 1 a 2 h ai b 2 K\ ^ ' S ' gai<i 2 ha 1 a 2 
a,i—la 2 ,b 2 ai — la 2 — l 

l,n/m l,m n/m l.m 

+Ka ^ ^ 9ai 02^6162 ~ Ko E E 9a ia2 h ai b 2 



ai,bi a 2 ,b 2 

If we take g = h 



a± = l a 2 ,b 2 



a<b 



tY; K ab\ga\ 2 = {Kx-K Q ) 

ai=] 

n 

-K x J2 \9a I 2 +^0 



12 = 1 

n 

E^ 



(A54) 



If we take g = h 
2^Kab\g a \ 2 



a<b 



n/m / m \ 2 

= (K X -if )E E W2 

ai=l \a 2 = l / 
n I n \ 2 

^E(f«) 2 +^ E?« ( A55 ) 



Substituting into_Eqs. (|A16IA17|) both Eq. (|A"54f - with 
K = X and g = h = - and Eq. (jA"55|) - with X = p, 
and g = h = e 1 ^ - one obtains 



rlRSB 



n/m m 

jj [] e-^M'i^] (A56) 

ai — 1 a2— 1 

/m rn 2 

^ e^-i^i (A57) 

a 2 = l 



AA x ^ 



ai — l 

Ai Ao 

n/rn i 



E' 



A/2 



EE' 

ai=l \a 2 = l 



a— 1 \a,= l 



+^E e2l0 + -E e 

a— 1 a— 1 

AA = Ai — Aq A/2 = ft\ — (2q 
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Using the identities 
5R[(a6) 2 ] = 5R[(a6) 2 ] = 
a A R | 3 | 2 /2 f d (R d 



(ab + abf 



(A58) 



2tt 



e -lcr/2 e «[C^ 9 ] (A59) 



S 2 ]/2 _ /"°° -x 2 /2„zSR[\/As] 



27T 



(A60) 



where A is complex and Ar is real, we can linearize the 
dependence on e 1 ^ in the partition function Eq. (1A56[) 
using Gaussian integral expressions: 



3? 



Ao 
2 



T 2-' 



If 



AA 



log / P[0]ex P 3? WA£- |polXy* 



Mo$>-^ 

o=l 



(A61) 



/ 



a 2 = l 



(A62) 



log / P[l]ex P 3? CiVAA-R - |A/z| e ^ ai 



a 2 = l 



.r 



m m 

* 2 = 1 £12 = 1 

where we defined the Gaussian measures: 
= -^e"^/ 2 



/2vr 



(A63) 

(A64) 
(A65) 



V\k] = Vp{^)Vp{C,{)Vp{x k ) 
Eq. (|A56I) becomes 

Z 1A5B = e -nAf /2 x 

J n { y p w y ^ n e£(0aia2; o,i) } 



(A66) 



ai — 1 



a 2 — 1 



£(V>;0,l) = H[e**(Ci^/AA«HA/Ij ( A 67) 

(-1)] 



cf. Eq. (|46j) . In the n — > limit the 'phase contribution' 
to the replicated free energy is: 

(A68) 



1 



lirn^log^ = f 
n->o n 2 

2?[0]log 1 2?[1] 



,£(^0,1) 



and the free energy is 

lim ( 

n->o 



/3$ = lim Girsb[Q sp ; A sp ] 

n— »0 
2^.2 r 



(A69) 



32 
-«i (ko| 4 



2 



l-(l-m)(| gi | 4 + H 4 ) 



\r\ 4 



(1 - m)(X 1 q 1 + /iin) + m(\ Q q + fx r ) 



m 



/if + //m 



lim — log Z 

7i->0 n ' 



1RSB 



Saddle point equations. Deriving G/n with respect to the 
parameters we obtain the twelve self-consistency equa- 
tions determining the order parameter values at given 
external pumping intensity and amount of disorder. 

• Deriving w.r.t. to Q's parameter we obtain the 
specification of Eqs. (|A18tiA21[) for each 1RSB 
replica matrix sector and for f and m, cf. Eqs. 

dUSOD, 



Ao,i — 


P 2 ° 2 , 
4 


Mo,i = 


4 


A = 


8 


v = 


P Jo , 
— 1 



9o,i 1 9o,i I 



l r o,i| r ,i 



Ifl 2 f 



(A70) 
(A71) 
(A72) 
(A73) 



• Deriving w.r.t ft and v we obtain Eqs. (|55p . where 
we define 



C£ = (cos 0) £ 



=>£(<£;0,i) 



s £ = (sinr/>) £ 



(A74) 
(A75) 



• Deriving G w.r.t. Ao.i and /xo,i and equating to zero 
we obtain Eqs. (|5HI54[) . after having integrated by 
part in the Gaussian measures. To help the non- 
expert reader to easily derive the self-consistency 
equations we exemplify the calculation of Eq. (|51l) . 



, dG 
! 9Af 



= 1 - (1 - m)qf 



(A76) 



- J V[0] (C?c £ + Cl S c) m l^\ R - |A M | 

The latter term can be simplified by integrating by 
part 



Vp{y)yF(y) 



Vp(y) 



dF{y) 
dy 



(A77) 



with y = Ci, Ci m E q- (]A77j) . yielding 

(Ci'cc + Clsc 
(cos 2 



AX R - |A/i| x (A78) 
(1 - m)c 2 c + sin 2 cj) - (1 - m)s|) 
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The self-consistency equation can thus be rewritten 
as, cf. Eq. (jijT]) . 

1 - (1 - m)qf = 1 - (1 - m) J V[0] (c 2 c + s 2 c ) m 
^ = «4>J + «4Uo " (A79) 



Eqs. (152II54[) ) are analogously derived. We notice 
that since from the equations dG/dX^ — one 
obtains l = the values of the q overlap are 
real- valued and so are the values of A. 



* Electronic address: luca.leuzzi@cnr.it, claudio.conti@romal. 

1 T. Schwartz, G. Bartal, S. Fishman, and M. Segev, Nature 
446, 52 (2007). 

2 G. Roati, C. D'Errico, L. Fallani, M. Fattori, C. Fort, 
M. Zaccanti, G. Modugno, M. Modugno, and M. Ingus- 
cio, Nature 453, 895 (2008). 

3 J. Billy, V. Josse, Z. Zuo, A. Bernard, B. Hambrecht, 
P. Lugan, D. Clement, L. Sanchez-Palencia, P. Bouyer, 
and A. Aspect, Nature 453, 891 (2008). 

4 C. Conti, L. Angelani, and G. Ruocco, Phys.Rev.A 75, 
033812 (2007). 

5 I. V. Shadrivov, K. Y. Bliokh, Y. P. Bliokh, V. Freilikher, 
and Y. S. Kivshar, Phys. Rev. Lett. 104, 123902 (2010). 

6 D. S. Wiersma, Nature Physics 4, 359 (2008). 

7 M. Bamba, S. Pigeon, and C. Ciuti, Phys. Rev. Lett. 104, 
213604 (2010). 

8 U. Bertolozzo, S. Residori, and P. Sebbah, ArXiv e-prints 
(2010), 1003.4931. 

9 T. Bienaime, S. Bux, E. Lucioni, P. W. Courteille, N. Pi- 
ovella, and R. Kaiser, Phys. Rev. Lett. 104, 183602 (2010). 

J. D. Bodyfelt, T. Kottos, and B. Shapiro, Phys. Rev. Lett. 
104, 164102 (2010). 

1 C. Conti, Phys. Rev. E 72, 066620 (2005). 

2 C. Conti, M. Peccianti, and G. Assanto, Opt. Lett. 31, 
2030 (2006). 

3 V. Folli and C. Conti, Phys. Rev. Lett. 104, 193901 (2010). 

4 S. E. Skipetrov, Phys. Rev. E 67, 016601 (2003). 

5 O. Zaitsev and L. Deych, arXiv:0906.3449 (2009). 

6 H. Cao, Waves in Random Media and Complex Media 13, 
Rl (2003). 

7 H. E. Tureci, L. Ge, S. Rotter, and A. D. Stone, Science 
320, 643 (2008). 

8 K. L. van der Molen, A. P. Mosk, and A. Lagendijk, Phys- 
ical Review A 74, 053808 (2006). 

9 R. G. S. El-Dardiry, A. P. Mosk, O. L. Muskens, and A. La- 
gendijk, Phys. Rev. A 81, 043830 (2010). 

C. Conti, M. Leonetti, A. Fratalocchi, L. Angelani, and 
G. Ruocco, Phys. Rev. Lett. 101, 143901 (2008). 

1 R. Weill, B. Levit, A. Bekker, O. Gat, and B. Fischer, Opt. 
Express 18, 16520 (2010). 

2 R. Weill, B. Fischer, and O. Gat, Phys. Rev. Lett. 104, 
173901 (2010). 

3 P. Suret, S. Randoux, H. R. Jauslin, and A. Picozzi, Phys. 
Rev. Lett. 104, 054101 (2010). 

4 U. Bortolozzo, J. Laurie, S. Nazarenko, and S. Residori, J. 
Opt. Soc. Am. B 26, 2280 (2009). 

5 S. K. Turitsyn, S. A. Babin, A. E. El-Taher, P. Harper, 
D. V. Churkin, S. I. Kablukov, J. D. Ania-Castanon, 
V. Karalekas, and E. V. Podivilov, Nat Photon 4, 231 
(2010). 

6 H. A. Haus, IEEE J. Quantum Electron. 6, 1173 (2000). 

7 L. Leuzzi, C. Conti, V. Folli, L. Angelani, and G. Ruocco, 
Phys. Rev. Lett. 102, 083901 (2009). 



mfn.T 

— 29" 



40 
41 



44 
48 



50 
51 
52 

53 
54 



ii. Kirkpatrick and B. Selman, Science 264, 1297 (1994). 
R. Monasson, R. Zecchina, S. Kirkpatrick, and et al, Na- 
ture 400, 133 (1999). 

M. Mezard, G. Parisi, and R. Zecchina, Science 297, 812 
(2002). 

M. Mezard and A. Montanari, Information, Physics, and 
Computation (Oxford University Press, 2009). 

C. Conti and A. Fratalocchi, Nat. Physics 4, 794 (2008). 
H. A. Haus, Waves and Fields in Optoelectronics 
(Prentice-Hall, Englewood Cliffs, N. J., 1984). 

K. Sakoda, Optical Properties of Photonic Crystals 
(Springer- Verlag, Berlin, 2001). 

G. Hackenbroich, C. Viviescas, B. Elattari, and F. Haake, 
Phys. Rev. Lett. 86, 5262 (2001). 

G. Hackenbroich, C. Viviescas, and F. Haake, Phys. Rev. 
A 68, 063805 (2003). 

L. Angelani, C. Conti, G. Ruocco, and F. Zamponi, Phys. 
Rev. Lett. 96, 065702 (2006). 

L. Angelani, C. Conti, G. Ruocco, and F. Zamponi, Phys. 
Rev. B 74, 104207 (2006). 

A. Yariv, Quantum Electronics (Saunders College, San 
Diego, 1991). 

W. E. Lamb, Phys. Rev. 134, A1429 (1964). 

M. Weigt, R. White, H. Szurmant, J. Hoch, and T. Hwa, 

PNAS 106, 67 (2009). 

T. Mora, A. Walczak, W. Bialek, C. Callan, and G. Cur- 
tis Jr., PNAS 107, 5405 (2010). 

P. Meystre and M. SargentHI, Elements of Quantum Op- 
tics (Springer, 1998). 

D. S. Wiersma and S. Cavalieri, Nature 414, 708 (2001). 
M. Leonetti and C. Conti, J. Opt. Soc. Am. 27, 1446 
(2010). 

F. Dalfovo, S. Giorgini, L. Pitaevskii, and S. Stringari, 
Rev. Mod. Phys. 71, 463 (1999). 

H. Stoof, Journ. of Low Temp. Physics 114, 11 (1999). 
R. A. Duine and H. T. C. Stoof, Phys. Rev. A 65, 013603 
(2001). 

R. A. Duine, B. W. A. Leurs, and H. T. C. Stoof, Phys. 

Rev. A 69, 053623 (2004). 

T. M. Lifshitz, Adv. Phys. 13, 483 (1964). 

L. Leuzzi and G. Parisi, J. Stat. Phys. 103, 679 (2001). 

L. Angelani, C. Conti, L. Prignano, G. Ruocco, and 

F. Zamponi, Phys. Rev. B 76, 064202 (2007). 

A. Gordon and B. Fischer, Opt. Comm. 223, 151 (2003). 

J.-L. Barrat, J.-N. Roux, and J. -P. Hansen, Chem. Phys. 

149, 197 (1990). 

J.-P. Hasen and S. Yip, Transp. Theory Sta. Phys. 24, 
1149 (1995). 

W. Kob and H. Andersen, Phys. Rev. Lett. 73, 1376 
(1994). 

W. Kob and H. Andersen, Phys. Rev. E 51, 4626 (1995). 
W. Kob and H. Andersen, Phys. Rev. E 52, 4134 (1995). 
F. Sciortino, W. Kob, and P. Tartaglia, Phys. Rev. Lett. 



20 



7M 



83, 3214 (1999). 

60 M. Mezard and G. Parisi, Phys. Rev. Lett. 82, 747 (1999). 

61 B. Coluzzi, G. Parisi, and P. Verrocchio, Phys. Rev. Lett. 

84, 306 (2000). 

62 E. Marinari, G. Parisi, and F. Ritort, J. Phys. A 27, 7615 
(1994). 

63 E. Marinari, G. Parisi, and F. Ritort, J. Phys. A 27, 7647 
(1994). 

64 L. Cugliandolo, J. Kurchan, and G. Parisi, Phys. Rev. Lett. 
74, 1012 (1995). 

65 D. Sherrington and S. Kirkpatrick, Phys. Rev. Lett. 35, 
1792 (1975). 

66 G. Parisi, Phys. Rev. Lett. 43, 1754 (1979). 

67 G. Parisi, J. Phys. A 13, L115 (1980). 

68 M. Mezard, G. Parisi, and M. A. Virasoro, Spin glass the- 
ory and beyond (World Scientific, Singapore, 1987). 

69 M. Mezard, G. Parisi, and M. Virasoro, Europhys. Lett. 1, 
77 (1986). 

70 F. Guerra, Comm. Math. Phys. 233, 1 (2003). 

71 M. Talagrand, Ann. Math. 163, 221 (2006). 

72 M. Mezard, G. Parisi, N. Sourlas, G. Toulose, and M. Vi- 
rasoro, Phys. Rev. Lett. 52, 1156 (1984). 
T. Castellani and A. Cavagna, J. Stat. Mech. P05012 
(2005). 

74 L. Leuzzi and T. M. Nieuwenhuizen, Thermodynamic of 
the glassy state (Taylor & Francis, 2007). 

75 E. Gardner, Nucl. Phys. B 257, 747 (1985). 

76 A. Crisanti and H.-J. Sommers, Zeit. Phys. B 87, 341 
(1992). 

77 D. Thouless, P. Anderson, and R. Palmer, Phil. Mag. 35, 
593 (1977). 

78 A. J. Bray and M. A. Moore, J. Phys. C: Solid State Phys. 
13, L469 (1980). 

79 A. Crisanti, L. Leuzzi, G. Parisi, and T. Rizzo, Phys. Rev. 
B 68, 174401 (2003). 

80 A. Crisanti, L. Leuzzi, and T. Rizzo, Eur. Phys. J. B 36, 
129 (2003). 

81 A. Annibale, A. Cavagna, I. Giardina, and G. Parisi, J. 
Phys. A 36, 10937 (2003). 

82 A. Crisanti, L. Leuzzi, G. Parisi, and T. Rizzo, Phys. Rev. 
B 70, 064423 (2004). 

83 A. Crisanti, L. Leuzzi, G. Parisi, and T. Rizzo, Phys. Rev. 



Lett. 92, 127203 (2004). 

84 T. Aspelmeier, A. Bray, and M. Moore, Phys. Rev. Lett. 
92, 087203 (2004). 

85 A. Crisanti, L. Leuzzi, and T. Rizzo, Phys. Rev. B 71, 
094202 (2005). 

86 M. Muller, L. Leuzzi, and A. Crisanti, Phys. Rev. B 74, 
134431 (2006). 

87 S. Edwards and P. Anderson, J. Phys. F 5, 965 (1975). 

88 A. Gordon and B. Fischer, Phys. Rev. Lett. 89, 103901 
(2002). 

89 P. Ehrenfest, Proc. Roy. Acad. Amsterdam 36, 154 (1980). 

90 If Jo = we are in the completely disordered case Rj = oo 
(also realizable by means of a finite Jo and o" 2 , = oo). In 
Ref. [13] V has been defined as ^Pk B T batil , simply amount- 
ing to an adimensional rescaling V — > Pv^JoAmth w.r.t. our 
model case. 

91 A factor of 8 has to be considered because of the over- 
counting of terms in the Hamiltonian of the model stud- 
ied in Ref. [37l ] with respect to Eq. I]13[l. This factor can 
be absorbed into the temperature yielding the pumping 
threshold Ppml = \/8/T . If we insert To ^ 0.717, i.e., the 
temperature at which the FM phase first appears in com- 
plete absence of disorder we obtain Vpml — 3.34. This also 
exactly corresponds to the spinodal value of V = 3.3412 
for Jo/a j — > oo in the present model. 

92 Since in our model the dynamic variables are continuous 
phases, the whole low T phase is consistently described by 
the 1RSB solution, unlike models with discrete variables 
such as the Ising p-spin modet^ where a further transition 
occurs at the so-called Gardner temperature. 

93 Technically speaking, this is due to the fact that all the 
terms of the free energy functional depending on two repli- 
cas observables have n — 1 or n — m factors in front, cf. 
App. [X] and in the n — > limit this factors change. 

94 The very existence of a Kauzmann temperature, also called 
the ideal glass transition temperature, in structural glasses 
is, actually, a matter of debate. 

95 The paramagnetic phase exists as metastable also at T < 
T s but the phase space is disconnected and the ergodicity 
is broken because of infinite barriers. 

96 The case at Jo = was explicitly considered in Ref. [37l ]. 



